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Получена формула, связывающая электростатическую компоненту дав- 
ления набухания с плотностью заряда и степенью дисперсности лиофиль- 
ной дисперсной системы. Формула справедлива для произвольных потенциа- 
лов поверхности и геометрической модели системы при условии справедли- 
вости гомогенного приближения (приближения постоянства потенциала 
двойного слоя внутри дисперсии), для этого условия получены количествен- 
ные неравенства. 


Электростатическая компонента расклинивающего давления в дис- 
персных системах изучается экспериментально и теоретически уже не- 
сколько десятков лет. Трудами Дерягина и его школы установлено, что 
это явление играет исключительно важную роль почти во всех коллоид- 
но-химических системах с водной дисперсионной средой и во многих 
других жидких дисперсионных средах [1]. Этот факт — общепризнанное 
достижение коллоидной науки, вошедшее в учебники [2—4]. Ознаком- 
ление с характером теоретических исследований в области электроста- 
тического расклинивающего давления показывает, что основное внима- 
ние уделяется той области параметров модели двойного слоя (ДС), 
которая соответствует классическим объектам коллоидной химии: ла- 
тексам, частицам окислов металлов (лиофобным коллоидам), и зна- 
чительно меньшее число работ посвящено моделированию ДС и вычис- 
лению расклинивающего давления для ситуаций, которые соответствуют 
лиофильным объектам (полиэлектролиты, биополимеры, синтетические 
полимерные мембраны, иониты). Сравнительная характеристика этих 
двух классов объектов (не претендующая на точность числовых значе- 
ний) приведена в табл. 1. 

Анализ литературы показывает, что даже немногие имеющиеся ра- 
боты по моделированию давления набухания лиофильных систем не свя- 
заны с основной массой работ по теории расклинивающего давления. 
Эти области предстают в настоящее время как бы независимыми, тогда 
как на самом деле физические основы этих теорий одни и те же, а ма- 
тематический аппарат различается лишь постольку, поскольку нам 
необходимо работать с различными наборами параметров одной и той 
же модели. Такой отрыв не лучшим образом сказывается на уровне тео- 
ретического понимания явления набухания лиофильных заряженных си- 
стем. В частности, даже тот факт, что широко используемое в теории 
мембранных явлений соотношение Доннана вытекает как частный слу- 
чай из решения уравнения Пуассона — Больцмана, сравнительно мало’ 
известен среди специалистов по мембранам. 

Модель для описания давления набухания. Давление набухания, как 
и расклинивающее давление (здесь и далее мы рассматриваем только 
электростатическую компоненту), есть результат отталкивания фикси- 
рованных зарядов дисперсной фазы, экранируемых ионами, находящи- 
мися в растворе. Поэтому в теоретическом отношении эти явления тож- 
дественны, и для описания давления набухания должен использоваться 
аппарат теории расклинивающего давления. Основное отличие объектов, 
в которых наблюдается давление набухания, от объектов, уже подробно 
описанных теоретически, состоит в том, что первые имеют намного мень- 
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Таблица 1 
Особенности лиофобных и лиофильных дисперсных систем, существенные 
для моделирования расклинивающего давления 
Е ЕН ВЕНЕ ЕЕ ОН ЕВА 


Система 
Свой ство 
лиофобная лиофильная 

ЕВЕ ИЕ СЕНО ЗНС 
Размер частиц 2000 А > ЗА 
Объемная доля Малая (две частицы) До 70% 
Относительная толщина двой- | Малая 1 и более 

ного слоя по сравнению с 

размером частицы 
Перекрытие ДС Существенно только в облас-| Может быть существенно во 

ти контакта всем объеме 

Взаимное расположение час- | Случайное Относительно упорядоченное 

тиц 

° Основное кинетическое Коагуляция Пептизация 
свойство 


Форма проявления поверх- Замедление или предотвра- | Давление набухания 
ностных сил отталкивания щение коагуляции 


ший характерный размер (большую степень дисперсности), чем вторые 
(см. табл. 1). Это количественное различие приводит к качественному 
отличию проявления сил отталкивания: если при моделировании лиофоб- 
ных объектов основную роль играет взаимодействие между двумя ча- 
стицами, то в лиофильных — коллективное взаимодействие между мно- 
гими «частицами» (будь то реальными частицами полиэлектролита или 
условно выделяемыми сегментами матрицы мембраны). Такое отличие 
обусловлено тем, что типичный лиофобный объект велик по сравнению 
с толщиной диффузного ДС 

а>1 (1) 


где а — характерный размер объекта (здесь и далее все линейные раз- 
меры измерены в единицах дебаевского радиуса экранирования)\, а ти- 
пичный лиофильный объект меньше этой толщины 


а<1 (2) 


Таким образом, задача вычисления давления набухания сводится к рас- 
чету строения ДС в системе с толстыми и сильно перекрытыми двойны- 
ми слоями соседних частиц. Работы этого плана почти без исключений 
основываются на предположении о малости потенциала поверхности ча- 
стиц 6 (здесь и далее все потенциалы измерены в единицах ЮТ/Ё) 


15151 (3) 
Исходное уравнение Пуассона — Больцмана 
Аф=зВф (4) 


где ф — потенциал ДС (мы ограничиваемся случаем 1-1-электролита) 
при этом превращается в линейное уравнение Дебая 


Аф=ф (5) 


которое и решается для той или иной геометрической модели дисперсной 
системы [5—8]. 


1 Особый случай составляют предельно мелкодисперсные лиофобные коллоиды, 
например красный золь золота, размер которых соизмерим с дебаевским радиусом. 
Для них область перекрытия двойных слоев при контакте может быть очень большой 
и в то же время основное явление для них — коагуляция, а не пептизация. Такое по- 
ведение не укладывается в предложенную классификацию; расчет расклинивающего 
давления в этом случае весьма сложен [1, $ 6, 7; 24]. 
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Особое направление в исследовании расклинивающего давления в 
концентрированных дисперсных системах представляет цикл работ по 
вычислительному моделированию уравнения состояния суспензий заря- 
женных сферических частиц [9—12]. Для наших целей, однако, резуль- 
таты этого моделирования не имеют значения, так как исходным момен- 
том моделирования непременно является задание в явном виде потен- 
циала взаимодействия между частицами (в цитированных работах 
задавался парный потенциал из теории ДЛФО); тем самым возможные 
коллективные взаимодействия полностью игнорируются. 

Рассматриваемая нами модель основывается на уравнении Пуассо- 
на — Больцмана (4), но для вычислений мы будем использовать не 
предположение о малости потенциала поверхности частиц (3), а пред- 
положение о сильном перекрытии ДС соседних частиц. Это приближение 
будем называть гомогенным (так как потенциал ДС почти постоянен 
внутри такой дисперсной системы), или доннановским, так как выпол- 
няется соотношение Доннана для связи заряда ДС с потенциалом по- 


верхности 
[6 $ (2 
-—в(- +И а. (6) 


2, 


где С — объемная плотность фиксированного заряда; С, — концентрация 
электролита в растворе, находящемся в равновесии с дисперсией. Вна- 
чале рассмотрим вопрос о давлении набухания в этом пределе, а затем 
вопрос о пределах применимости самого приближения. 

Давление набухания в доннановском приближении. При вычислении 
давления набухания Р„ воспользуемся формулой Ленгмюра 


Р„=2ЮТС, (СН фв — 1) (7) 


где фз — потенциал ДС в точке системы, где электрическое поле ДС рав- 
но нулю. Используя (7), предполагаем, что условия, при которых фор- 
мула Ленгмюра дает верное выражение для расклинивающего давления, 
выполнены (подробно этот вопрос рассмотрен в [1, $ 6.1]). 

В предположении очень сильного перекрытия ДС, точный критерий 
которого будет установлен ниже, потенциал ДС постоянен всюду внутри 
дисперсной системы 

ф(Г) =6 


Фв=б (8) 


и таким образом, расклинивающее. давление очень просто выражается 
через 5-потенциал такой системы (называемой также доннановским по- 
тенциалом) 


Р.=28ТС, (св — 1) (9) 


Важность формулы (9) состоит в том, что она дает выражение для дав- 
ления набухания любой дисперсной системы независимо от ее геометри- 
ческого строения: две находящиеся в контакте дисперсные системы с 
одинаковыми доннановскими потенциалами будут иметь одинаковое дав- 
ление набухания при условии, что для обеих систем справедливо гомо- 
генное приближение. Когда же две дисперсные системы будут иметь оди- 
наковый доннановский потенциал? Чтобы ответить на этот вопрос, вы- 
разим б-потенциал через плотность заряда и геометрические параметры. 
Так как потенциал ДС постоянен, то плотность объемного заряда в нем 
также постоянна и равна р 


р=—2ЕС, 366 (10) 
Общий заряд диффузной обкладки ДС равен @ 
9=ру=—2РС.У 3,5 (11) 
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где У — объем порового пространства системы (так как в гомогенном 
пределе ДС занимает весь предоставленный ему объем). Этот заряд. 
уравновешивается фиксированным зарядом дисперсной фазы, равным 


—0@=00$ (12) 
где с — поверхностная плотность заряда, $ — площадь поверхности дис- 
персной фазы. Приравнивая (11) и (12), получаем 


0$ 
Су (13} 


Таким образом, 5-потенциалы систем будут равны, если у них равны па- 
раметры безразмерного заряда 


(27) = ( т) — Чет (14) 
СУ 1 СУ 2 


Чтобы практически воспользоваться критерием (14), предположим, что. 
концентрации электролита одинаковы 


С, =14ет (15) 


и одинаковы плотности заряда 
о=1!4ет (16} 


Тогда равенство 5-потенциалов будет достигнуто при условии равенства 
геометрического параметра 


А=-И = 14ет (17) 


независимо от конкретной формы поверхности дисперсной фазы. 

В чем смысл параметра ^? Он тесно связан с таким популярным по- 
казателем дисперсности, как удельная поверхность. В самом деле, обо- 
значим через 5,, площадь поверхности дисперсной фазы на единицу 
веса сухой, дисперсной фазы. Если плотность дисперсной фазы обозна- 
чить через р;, а объемную долю ее — через «, то 

5 оби 1 
+ же С ит к о ет 


Уорда У ра № @ 
откуда 


ды | (18) 


(У, — общий объем дисперсной системы). 

Таким образом, показатель дисперсности /, тесно’ связан с традицион- 
ным показателем $,„, но связь эта не непосредственная, так как объем- 
ная доля и, возможно, удельная поверхность лиофильной системы изме- 
няются в процессе набухания. Поэтому (18) представляет собой не 
столько расчетную формулу для ^, сколько иллюстрацию взаимосвязи 
величин Хи $,д; точным же определением параметра является (17) 


2, объем дисперсионной среды 
— —6 о 
площадь поверхности дисперсной фазы 


Для теоретического исследования можем задаться определенной гео- 
метрической моделью дисперсной системы и, вычислив для нее это отно- 
шение, устанавливать эквивалентность” различных геометрических мо- 
делей. Некоторые результаты такого рода представлены в табл. 2. Из 
нее, например, видно, что результат моделирования дисперсной систе- 
мы сферическими полостями будет совпадать с результатом моделирова- 
ния цилиндрическими капиллярами при условии, что радиус капилляров 
составляет */з радиуса полостей и т. д. 
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Таблица 2 


Свойства различных геометрических моделей 


Исходный гед- | Объем дис- Ра- Площадь Геометри- Решение крае- нь а 
Свойство метрический персионной диус | поверхности ческий Точка В Точка А вой задачи потенц, Функция 
параметр, а среды га ны ре в (20) в (19) ее к (23) Фи(В) ры [ (@) в (24) 
Сферическая по-| Радиус лав — Ала? а Центр по- | Поверх- 4-5 я 21. 3 
лость 3 . 3 лости ность дис- | в (22 — а?) а" 
: персной 
фазы 
Цилиндрический > ла?[, —. 2лаГ, а Центр ка- |То же я . 1 И 
капилляр га пилляра. © — а’) та 
Плоский капил- | Полуширина 2аЁН — 2ЕН а Центр ка- » й р 1 1 
ляр | | пилляра Ру (* — а?) НЕ = 
Решетка цилинд-| Радиус ци- |л (53 —а)Ё| а 2лаЁ а ({— о) | Поверхность › .р8 Е т а? 
ров линдра Иа 2 а | ячейки д 4 + Ша) 
5 ша’ 
= Ид а 1-2 ша 
2 2% [е2 
а? 
—= + 
р2 
- РЗ ша 
Решетка сфер | Радиус сфе- | 4л 8) а 4? |а (1—а)| То же » й. фи > а 9 (1 
ры 3 ( а $ а 6 35 20; (1 — о) \ 6 
| м 4 1 1 1 
6 34а 6 3 а Ва ) 


ы 5 


Рис. 1. Системы координат и исходные геометрические параметры 
моделей: 1 — сферическая полость; 2 — цилиндрический капилляр; 
3 — плоский капилляр; 4 — решетка цилиндров; 5 — решетка 


сфер 


Таким образом, все геометрические модели с одинаковым парамет- 
ром А равноценны в том отношении, что предсказывают одинаковое дав- 
ление набухания. Этот результат важен, так как обосновывает право- 
мерность применения регулярных геометрических образов для модели- 
рования иррегулярного порового пространства лиофильной дисперсной 
системы. 

Пределы применимости доннановского приближения. Чтобы опреде- 
лить пределы применимости гомогенного описания дисперсной системы, 
необходимо решить уравнение Пуассона—Больцмана (4) с краевыми 


условиями 
Ф[4=5 (19) 
д 
| = (20) 
В 


где А —точки, лежащие на поверхности дисперсной фазы, В — точки, 
где по соображениям симметрии электрическое поле должно отсутство- 
вать. Рассмотрим пять вариантов геометрии (см. табл. 2): сферические: 
полости, цилиндрические капилляры, плоские капилляры, решетку ци- 
линдрических нитей и решетку сферических частиц (рис. 1). Этот набор: 
моделей охватывает практически все возможные регулярные геометри- 
ческие структуры. Модели, включающие объемы эллипсоидальной фор- 
мы, при некоторых предельных значениях параметров вырождаются в. 
рассматриваемое нами множество моделей (табл. 3), поэтому можно ду- 
мать, что их рассмотрение, весьма громоздкое в математическом отноше- 
нии, не принесет много нового. _ 

Из решения краевой задачи (4), (19), (20) мы должны найти поправ- 
ку к приближенному выражению для потенциала в точке В [8] 


9в =6- 9: (В) (21) 

Предполагая поправку ф. малой величиной, разложим уравнение (4) 
по ф: 

Дф: = 15 - О ($1) (22} 
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Таблица $ 


Предельные случаи геометрических моделей с эллипсоидальными поверхностями 


р рррХрХЁАЁ&АяЁ&АЁа©&р ЖЖ ЖЖ ЖЖЖЖЖЖЖЖЖпррЖЖЖЖЖЖЖЖЖЖЖпЖЖЖпррррЖЖЖЖЖЖЖЖЖЖЖЖЖЖпаяпяждпХао Христиан, 


Эксцентриситет эллипсоида вращения 


Модель 
стремится стремится к бесконечности 

Эллипсоиды-частицы: 

вытянутые Решетка сфер Решетка цилиндров 

сплюснутые То же Плоский капилляр 
Эллипсоиды-полости: 

вытянутые Сферическая полость Цилиндрический капилляр 

сплюснутые То же Плоский капилляр 


для которого краевое условие (19) заменится на 
ф:/д=0 (23) 


Решения краевой задачи (23), (22), (20) приведены в табл. 2. Эти ре- 
шения можно использовать для получения областей применимости го- 
могенного описания. Задавшись относительной погрешностью давления 
набухания = (в нашем примере г=5%), мы определяем эти области не- 


равенством 
1 = 
_ Из (@) ы: 
св г 


где выражения для {[(*) для разных геометрий приведены в табл. 2. 
Границы областей в зависимости от объемной доли приведены на рис. 2, 


Рис. 2. Области справедливости го- 
могенного описания концентрирован- 
ных дисперсных систем (они лежат 
ниже соответствующих кривых): 1-— би 
5—6=2; Г—5’—6=5; цифры у 
кривых отвечают различным геомет- 
рическим моделям в соответствии с 
нумерацией рис. 1 
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из которого видно, что если неравенство (24) выполняется для сфериче- 
ских полостей, то оно практически всегда выполняется и для всех других 
геометрий. Поэтому в качестве критерия справедливости гомогенного 
описания можно использовать неравенство 


1 = 
< (25) 
в ь 


Сравнительно жесткое ограничение на величину геометрического пара- 
метра физически обусловлено тем, что гомогенное описание справедливо 
в условиях перекрытия внутренних областей диффузной обкладки ДС, 
т. е. областей, где происходит быстрое снижение потенциала до величи- 
ны порядка единицы. Эта область существенно тоньше дебаевского ра- 
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диуса и с ростом 5-потенциала все больше утоньшается, в связи с чем 
неравенства (24), (25) содержат экспоненту &-потенциала. 

Очень близкие результаты получены в [13], правда только для ре- 
шетки сферических частиц. Анализ влияния полидисперсности формы 
частиц на расклинивающее давление, насколько нам известно, не про- 
водился. Работы по решению нелинейного уравнения Пуассона—Больц- 
мана в концентрированных дисперсных системах [14—19] по большей 
части не затрагивают вопроса о расклинивающем давлении, а исполь- 
зование получаемых авторами решений затруднено различными обстоя- 
тельствами математического характера. 

Особо следует остановиться на цикле работ [20—23], в котором элек- 
‘тростатическое расклинивающее давление в концентрированной дисперс- 
ной системе исследуется именно в рамках гомогенного приближения. 
Однако они, к сожалению, изобилуют физическими несообразностями и 
математическими ошибками. В качестве примера физически неверного 
зывода приведем окончательный результат двух статей [21, 22], в ко- 
торых утверждается, что одноименно заряженные частицы все же при- 
тягиваются (с оговоркой, что это относится только к концентрированным 
дисперсным системам, а две изолированные частицы с теми же поверх- 
ностными свойствами — отталкиваются). Вопросами о том, при какой 
объемной доле и по какой причине происходит столь удивительное из- 
менение свойств электрических сил, авторы не задаются. В качестве при- 
мера математической софистики укажем, что краевое условие (20) зада- 
ется на поверхности ячейки, содержащей две дисперсные частицы. Точ- 
ная форма этой поверхности неизвестна, навязывание ей простой 
геометрической формы менее привлекательно, чем для ячейки с одной 
частицей. Авторы не определяют явно поверхности, на которой они за- 
дают (20), и поэтому не имеют необходимого числа уравнений для ко- 
эффициентов ряда Фурье. Для обхода этой трудности они предполагают 
справедливыми некоторые нетривиальные соотношения между искомы- 
ми коэффициентами, полагая, что дифференцированием полученного 
явно решения можно определить форму поверхности для граничного 
условия (20). В результате таких ошибок большое количество громозд- 
ких вычислений, которое можно было бы использовать для определения 
пределов применимости гомогенного подхода, оказалось выполненным 
напрасно. 
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